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Abstract 

Scalar Tensor Vector Gravity(STVG) is one of modified gravity theories developed by 
John Moffat(2005). MOG is abbreviated name for this theory. This theory is added 
a massive vector field to Brans-Dicke theory. It can explain a galactic rotation curve 
and the structure formation without dark matter. It can also explain acceleration 
universe without dark energy. However,these are claims by the developer and 
collaborators. This theory was only inspected by simple approximate calculation. 
Therefore it needs more objective verifications. We will carried out verification 
from the viewpoint of N-body simulation. Such study is already accomplished 
by Brandao(2010). However, they did not precisely formulate N-body simulation 
on MOG. This paper shows formulation of the N-body simulation on MOG more 
precisely. 

^ ■ 1 Introduction 

The standard model of cosmology today, the ACDM model, provides an excellent fit to cosmological 
observations. 

But, most of the composition of the universe is invisible and undetectable which is named dark energy 
or dark matter. This fact provides a strong incentive to seek alternative explanations that can explain 
cosmological observations without dark matter or dark energy. 

Modified Gravity (MOG) [J ] has been used successfully to account for galaxy cluster masses [8], the 
rotation curves of galaxies [5] [6], velocity dispersions of satellite galaxies, and globular clusters [11]. It can 
explain the observation of the Bullet Cluster [G] without cold dark matter. Besides, MOG also meets the 
challenge posed by cosmological observations. In the paper(arXiv:0710.0364)[4],it is demonstrated that 
MOG produces an acoustic power spectrum, a galaxy matter power spectrum, and a luminosity-distance 
relationship that are in good agreement with observations. 

However, these are claims by the developer and collaborators. Therefore it needs objective verifications. 
We want to carry out verification of MOG from the viewpoint of N-body simulation. But, such study 
is already accomplished by Brandao et al. (2010) [7]. The result was negative for the Moffatian model. 
They claim that the Moffatian potential cannot maintain exponential disks in dynamical equilibrium 
and therefore cannot be consistent with observations. However, they did not precisely formulate N-body 
simulation on MOG. This paper shows formulation of the N-body simulation on MOG in comformity 
with Moffat et al. (2009) [13]. 



2 Theory of MOG 
2.1 Action 

The action of Moffat theory is constructed as follows [1]. STVG is formulated using the action principle. 
In the following discussion, a metric signature of [+, — , — , — ] will be used; the speed of light is set to c = 1, 
and wc arc using the following definition for the Ricci tensor: R^ = d a T« v - d v T™ a + T^T^ - T^T^. 
We denotes the Einstein-Hilbert Lagrangiant^G' = ~ 16 * G (i? + 2A) y/—g, where R is the trace of the 
Ricci tensor, G is the gravitational constant, g is the determinant of the metric tensor g^ u , while A is 
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the cosmological constant. Introducing the Proca action Maxwell-Proca Lagrangian for the STVG vector 
field <h„: 
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where = c^t/v — d u (f>^, (i is the mass of the vector field, uj determines the strength of the coupling 
between the fifth force and matter, and is a self-interaction potential. The three constants of the 
theory, G, /i and u, are promoted to scalar fields by introducing associated kinetic and potential terms 
in the Lagrangian density: 
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where V M denotes covariant differentiation with respect to the metric g^, while Vq, V^, and V u are the 
self-interaction potentials associated with the scalar fields. The STVG action integral takes the form 
S = J (Cg + £<p + Cs + Cm) d 4 x, where Cm is the ordinary matter Lagrangian density. 



2.2 Weak field approximation of MOG 

The field equations of STVG can be developed from the action integral using the variational principle. 
First a test particle Lagrangian is postulated in the formXxp = —m + autq^^^u^ , where m is the test 
particle mass, a is a factor representing the nonlinearity of the theory, q§ is the test particle's fifth-force 
charge, and it M = dx^/ds is its four-velocity. Assuming that the fifth-force charge is proportional to 
mass, 95 = Km, the value of n = WGnJui is determined and the following equation of motion is obtained 
in the spherically symmetric, static and weak gravitational field of a point mass of mass M: 
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l + a-a(l+r/A)e- r / A l , (3) 



where Gjv is Newton's constant of gravitation. Further study of the field equations allows a determination 
of a and A for a point gravitational source of mass M in the form: A = 1/fj, = ^jf-,a = The 

constants D and E various astronomical observation yield the following values: D ~ 6250M Q ' kpc , 
E ~ 25000Mq /2 , 

In the weak-field approximation, STVG produces a Yukawa-like modification of the gravitational 
force due to a point source. Nearby a source gravitythe repulsive force which comes from vector field 
counteracts attraction. But far from a source gravity, Yukawa-like repulsive force cannot reach because 
vector field is massive. 

Intuitively, this result can be described as follows: far from a source gravity is stronger than the 
Newtonian prediction, but at shorter distances, gravity is comparable to Newtonian. 

MOG advocate that we recognize this strong gravity appearing by leaving for far away as dark matter. 
[1 + a — a(l + r/A)e _r / A ] in Eq.3 is able to interpreted as effective gravitational constant. Fig.l shows 
relation between mass scale / distance scale and effective G. The maximum of effective G and the effective 
range grow large depending on the mass of object. (See Fig. 2) Therefore it can explain an effect of the 
dark matter of every scale. 

The reason why the MOG can explain galactic flat rotary curves without dark matter halo is as 
follows. Fig. 3 shows effective gravitatonal constant of lO lo M0(It is comparable to galaxy disk). On 
galaxy disk scale ,G e // is in proportion to a radius. Let us remember equation of Kepler motion. The 

revolution velocity of object having the circular orbitis is discribe as v = \J ^r-- If Geff is in proportion 
to r, velocity is constant. 

The above is qualitative explanation, but agreement of observed rotation velocity by Doppler effect 
is identified as predicted rotation velocity by brightness in 100 galaxy [ r >] [G]. 
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Figure 1: Relation between scale and effective G. X-axis shows a distance scale. Y-axis shows mass of the 
gravity source. 



Object Mass Varying effective G MaxG=(l+ct) 




Figure 2: Effective G and distance from the object. The example of three different weight's objects: 
1O 9 M O ,1O 1O M ,1O 13 M O . 
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3 Techniques and detail of N-body simulations 

The gravity acceleration equation of the STVG in a weak gravitational field is shown in the Eq.3. But, This 
cqutaion is adaped only gravitational field of "a point mass" . It is not applied simple superposition princi- 
ple of forces without self-contradiction. The precedent study by Brandao [7] treat MOG as simple Yukawa 
gravity. Yukawa-effective range is fixed as 13.96fcpc. But, MOG is not simple Yukawa gravity. Yukawa- 
effective range(=A)and increase of effective G(=a) depends on neighboring mass distribution. Every 
combination of an interacting particle(i-j) have Xji and ctji. 

The mutual gravity between the particle by STVG is written in Moffat et al. (2009) [13]. I let an 
equation written in the paper become disintegration. And the following expressions are derived by 
adapting itself to many body system. 

EG N mj(fl-fj) 
7i FT — k^ Geffihj) (4) 
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where £ is a parameter which is decided by observation. Even if the reader watches a numerical formula, 
the image of the physical meaning is hard to appear. In the next page, a specific example is explained 
using a figure. 
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Figure 4: Example of the acceleration calculation 

We think about the placement of the particle such as the Fig. 4. the gravitation that i particle gives 
to j particle is described in the following expressions: 



G N mj(fl-fj) 

a iJ = -j, h~r — k^ G eff(hj) 
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The calculation of gravitational acceleration on Moffat gravity needs the information of all the particles 
whenever we calculate force between two particles. That is why the calculation number of times becomes 
C(A r3 ).By the normal method, we cannot carry out such an enormous calculations. 
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We show below the method of reducing calculation number of times. We pay attention that M e ff 
which decide the constant of gravitation increment(= a) and Yukawa- effective range(= A) depends on 
"absolute value" of the distance between the particle. M e tf integrated a function of the mass density 
per unit of area of the shell around the gravity source. We can reduce calculation number of times by 
defining to an individual particle with the "mass shell" which is a state in Fig. 5. In the following, we 
abbreviate "mass shell" to "M shell" . 




Mshell[l]= - ffi 

Mshell[2]= tj|^-_- ffi 
Mshell[3]= jij^ Qsj js 

Mshell[4]= ^Olpc^^l 



Figure 5: Image of mass shell 

G e ff don't depend on direction from the gravity source particle to the gravity receive particle ,but 
only distance from source particle to receive particle. Therefore, we can define the G e //'s shell(G shell) 
concentric from the gravity source particle. We propose the following numerations from these hints. In 
the case of normal N-body simulation, we calculate vector sum of net force on every gravity receive paritclc 
from all the other gravity source particles. But, in the case of MOG,we calculate in the procedure that 
every gravity source particle gives force for all the other gravity receive particles. 

1. Getting mutual coordinates of all the other gravity receive particles from the gravity source parti- 
cle. In parallel, defining "mass shell". (The calculation number of times is N — 1) 

2. Defining "G shell" by using "mass shell". It is necessary sum of all the "mass shell" to calculate one" 
G shell" (The calculation number of times is Nmassshell x NGshell) 



3. The gravity source particle giving force for all the other gravity receive particles. Using "G shell" 
for this calculation. (The calculation number of times is N — 1) 



We carry out this calculation for all the gravity source particles repeatedly. (The calculation number 
of times is N) 
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Assuming that both "mass shell" and "G shell" have equal interval shell, the basic equation is renewed 
below. We define SM shell as interval of "mass shell" , and SGshell@&s interval of "G shell" . 
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G e f f shell (i, nGshell) = 
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where LXJ means rounding off fractions of X.For example, [1.2 J = 1, |_3.6j = 3. In reality, "mass shell" 
and "G shell" does not necessarily have to be equal interval. Considering function form of M e //,"mass 
shell" of interval should be small in a center area. Because the number of shell which wc discribc in a 
program code is not infinity ,this method can adopt for only the system whose size does not change. We 
treat a system of the dynamics equilibrium in most studies. Thus,we do not have any problem in this 
matter. 



a(j3)+=a newton (j3) x Gshell[4] rg hen ^ 




a(jl)+=a newton (jl)xGshell[2] 



a(j2)+=a newton (j2)xGshell[3 



Figure 6: Image of the acceleration calculation using "G shell' 
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The precision of this approximate calculation depends on the number of "mass shell" and "G shell" . 
To examine dependence of the precision by number of "shell" ,we calculated the acceleration which each 
particle of the cxpncntial disk received. The following is parameter of the disk.Thease are comparable to 
the Milky way. 

lR(radius) 25kpc 
lM(mass) 5 x 1O 1O M 
lT(time) 0.26Gry 
lV(vclocity) 92km/s 

Table 1: The standardization of the unit 



Rd (scale length of disk) 


2.5kpc 


R C ut (radius of disk) 


25.0kpc 


Zd (scale length of disk height) 


0.5kpc 


Ma (total mass of disk) 


5 x 10 10 M Q 


cps (softening length) 


0.048kpc 


Nbody(number of particles) 


10000 


£(MOG parameter) 


12.0 



Table 2: Setting parameters 



The radius of most outer "shell" (:5Mshell x NMsheii = SGshell x Ncsheii) is about 3.0. This length 
is enough to caluculatc acceleration of the pair of particle which have the longest distance in the disk. 
We examine the 7 model: N MsheU x N GsheU = 10 2 , 50 2 , 100 2 , 250 2 , 500 2 , 1000 2 , 2000 2 . For a comparison, 
we calculated the model of the 0(N 3 ) exact calculation. 

Then, we calculated the MOG acceleration of each particle in this disk. The result is shown in Fig. 7. 
This figure shows the ratio of central acceleration by the modified gravity and the central acceleration 
by the Newton gravity as an effective gravitational constant. The data plots in the Fig7 overlaps ,we 
cannot distinguish each model. The Fig8 expresses the lines which averaged each point in Fig. 7. Each 
line in the Fig8 also overlaps ,we cannot distinguish each model. However, it's the important that all 
other line except "10 2 shell model" overlaps with exact calculation model which has calculation number 
of times 0(N 3 ). After all, we don't need to divide interval of "shell" so much. If Nmassshell x NGsheU 
is smaller than iV,thc total calculation number of times becomes 0(N 2 ). In this case, the calculation 
load become equivalent to the normal N-body calculation. We are groping for the method to reduce the 
calculation number of times more. 

4 Postscript 

This paper is an interim report. Now,we are carrying out a large-scale calculation about galactic dynamics 
and large scale structure of universe. Nurmcrical calculation are carried out on Cray XT4 at Center for 
Computational Astrophysics(CfCA) of National Astronomical Observatory of Japan and SR16000 at 
YITP in Kyoto University. If a meaningful result of the simulation appears, we intend to submit an 
official article. 
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Figure 7: The ratio of central acceleration by the modified gravity and the central acceleration by the 
Newton gravity. 
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Figure 8: lines which averaged each points in Fig 7. 
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